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ABSTRACT 



The goal of this thesis is to investigate the work output of 
a heat driven thermoacoustic prime mover above onset of self- 
oscillation. The exponentially growing sound wave, generated when 
a prime mover is initially "turned on," was digitally sampled for 
a helium filled prime mover at pressures ranging from 238 kPa to 
500 kPa and at temperature differences ranging from onset to 400 K. 
This data was then digitally filtered by a 100 Hz band pass filter 
centered on the prime mover's fundamental frequency. A least mean 
squares fit was applied to the envelope of the filter's output in 
order to determine the temporal absorption coefficient ]3. From 13 , 
the quality factor was computed. These quality factors were then 
compared with thermoacoustic theory. The agreement between the 
theoretical predictions and the measured results is extremely good 
at high mean gas pressures. As the mean gas pressure decreases, 
however, the agreement between the theoretical value of onset and 
the predicted slope of the data, increasingly deviate. 
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INTRODUCTION 



As thermoacoust ics is a relatively new technology, its 
potential uses are still being investigated. One such use, 
which this thesis intends to precede, is the use of a 
thermoacoustic prime mover as a low frequency transducer. 
This possibility has been discussed by S. L. Garrett and T. B. 
Gabrielson as a viable alternative to the conventional sound 
sources used today. 

For conventional sources, the motion of a solid surface is 
required to generate a volumetric velocity. The pressure 
amplitude is consequently limited by the maximum displacement 
of the solid. The relationship between the pressure amplitude 
P and the angular frequency of oscillation w for a spherical 
source in an infinite medium is 

P = Q-^- = Q-^-> U) 

2 kr 4 nr 

where the volumetric velocity Q is given by 

Q - f vndA. (2) 

J surface 



Also, p is the density of the medium, a is the speed of sound 
in the medium, X is the wavelength of the oscillating wave, r 
is the range from the source, and v, n, and A are the surface 
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velocity, the unit vector normal to the surface, and the 
surface area of the source, respectively. Assuming sinusoidal 
oscillations, the surface velocity can be expressed as v = cod, 
where d is the displacement of the surface. Consequently Q is 
proportional to co, the pressure amplitude to cj 2 , and the 
radiated acoustic power to co 4 . As a result, as frequency is 
decreased, the surface area must be increased quadratically 
in order to maintain a constant pressure amplitude in the far 
field (P oc co 2 *A) . Thermoacoustics offers a solution to this 
problem because the driving force can be applied directly to 
the fluid instead of through an oscillating surface. With 
only fluid motion, the limitation to the velocity no longer 
comes from the elastic constraints of a vibrating solid but 
from the hydrodynamic (e.g., turbulence) limits of the fluid. 
(Garrett et al., 1990, pp. 162-163) The actual phenomenon of 
thermoacoustics has been discussed at length elsewhere. G. W. 
Swift's 1988 article "Thermoacoustic engines", for example, 
gives a very thorough coverage of the subject. It therefore 
will not be explained in this thesis. 

Thermoacoustics research at the Naval Postgraduate School 
has been directed at thoroughly verifying thermoacoustic 
theory. In 1989, Hsiao-Tseng Lin studied a thermoacoustic 
prime mover below onset (Lin, 1989) . The conclusions of his 
thesis indicated that theory based on short-stack and boundary 
layer approximations did not hold valid in all cases. He also 
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found that above onset, the acoustic amplitude quickly became 
quite large, and the waveforms highly nonlinear. 

The next logical progression in this research is to study 
a prime mover above onset while the acoustic amplitude is 
small and to apply a non short stack, boundary layer theory. 
It is intended that the information gained will help lead to 
an understanding of prime movers at large amplitudes. The 
purpose of this thesis is to take the first step by looking at 
the initial growth of the oscillations as they build up to 
steady state. 
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II 



THEORY 



A. INTRODUCTION 

A prime mover can be divided into five major regions. 
These five regions, as shown in Figure 1, are the ambient end, 
the ambient heat exchanger, the stack, the hot heat exchanger, 
and the hot end. In typical configurations, the ambient and 
hot ends are empty capped tubes maintained at constant 
temperature. Since we chose to apply a differential 
temperature to the prime mover by heating one end while 
keeping the other end at ambient temperature, hot and ambient 
are used to describe the different regions. The ambient and 
hot heat exchangers both consist of a series of plates used to 
maintain the opposing ends of the stack region at their 
respective temperatures. Finally, the stack region is where 
the conversion of thermal to acoustic energy takes place. 

The purpose of this thesis is to compare the measured work 
output of a prime mover above onset of self-oscillation to the 
theoretical value. Unfortunately, it is not a simple matter 
to directly measure the work output. What is easily measured 
though, is the acoustic signal generated by the prime mover. 
The evolution of this signal can be related to the work output 
as discussed below. 
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When a resonator is excited and then allowed to decay, the 
amplitude of the oscillations decay according to 

g-ut/2 

where t is time, Q is the quality factor of the resonator and 
j 2 = -1. (Swift, 1988, p. 1162) In terms of the temporal 
absorption coefficient j 8, this exponential behavior is 

e -pt e jut f (4) 

where 0 = w/2Q. The initial buildup of oscillations in a 
prime mover can be characterized by a negative value of (S . 
Thus, the amplitude initially grows exponentially. As the 
oscillations become larger in amplitude, however, the value of 
the net absorption coefficient changes and eventually equals 
zero. This results in steady state oscillations of finite 
amplitude as can be seen in Figure 2. It is the initial 
growth of the signal that is of interest to us. 

) 3 , which is equivalent to the inverse of the rise time, 
can be directly obtained from analyzing the envelope of the 
acoustic signal generated by the prime mover. Along with 13, 
the frequency of oscillation is also measured so that 1/Q may 
be computed. Recall that 1/Q = 2/3/w = (3/nf, where f is the 
oscillation frequency measured in hertz. The choice of using 
1/Q as the parameter of interest comes from the fact that it 
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Signal sampled at O.Olsec/pt (For example only) 

0 2 1 1 , 1 i , 




- 0 . 2 1 * * 4 ' * * * 

0 5 10 15 20 25 30 35 40 



Time (sec) 



Figure 2 Signal which shows the changing of 0 from an 
initial negative value (evidenced by the exponential 
buildup) , to zero (as indicated by the steady state 
result) . 
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is directly proportional to the work output of the prime 
mover. As discussed by Swift, 1/Q relates to the average 
dissipated acoustic power E of a prime mover by 



1 _ E 

Q 0 )E st ' 

* 



(5) 



where E ST is the energy stored in the resonator (Swift, 1988, 
p. 1162). Rather than using the dissipated acoustic power E, 
we prefer to work with the generated acoustic power W. W is 
commonly referred to as the work output. By convention, W is 
positive when power is generated while E is positive when 
power is dissipated. Therefore, E = -W, and 



1 _ — Pv^ 

Q a >E st ' 



( 6 ) 



B. AVERAGE GENERATED ACOUSTIC POWER/ W sUck 

The acoustic power generated in the stack section of a 
prime mover cannot flow or escape in a direction perpendicular 
to the stack plates, the y-direction. Energy transport is 
allowed however, in the x-direction which is parallel to both 
the plates and the temperature gradient. A measurement of the 
difference in average acoustic intensity between the two ends 
of the stack can therefore be used to determine the acoustic 
power generated by the stack. 
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The acoustic power generated by a short segment of stack 



is given by 



w. 



stack 



S [ ( pu ) 



right 
end of 
segment 



(pu) left. ] 
end of 
segment 



= Sax-^- (pu) 

dx 



(7) 



where S is the cross sectional area of the fluid in the stack, 
pu is the time average acoustic intensity, and Ax is the 
length of the short segment. Note that SAx is the total 
volume of the fluid in the stack region. By dividing through 
by this volume, the power per unit volume is obtained; 



W. 



w = 



stack _ 

Sax 



d , — \ 1 _ r d {Qj} 

-(pu) - -Re [ Pl — 



dp l ! - 
dx 



( 8 ) 



where a tilde denotes complex conjugate, and p, and u, 
represent the spatial dependence of the pressure and velocity, 
respectively, (i.e. p(x,t) = p^" 1 and u(x,t) = u,^) . Also, 
the angled brackets denote the spatial average of velocity 
across the fluid gap between the plates. It is necessary to 
spatially average the acoustic velocity across the fluid gap 
because, with the presence of viscosity, u is also y- 
dependent. Expressions for the four terms on the right hand 
side of this equation will be derived next in order to arrive 
at a complete equation for work output. (Swift, 1988, p. 1159) 
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Since the prime mover is essentially a pipe with rigid 
ends, the pressure behaves sinusoidally like a cosine wave. 



= P A cos(kx) 



(9) 



where P A is the peak pressure amplitude, k is the wave number, 
and x = 0 at the capped end of the ambient end. 

Similarly, it can be shown that the velocity in the stack 
is represented by 



where p m is the mean density of the fluid, 1 is one-half the 
plate thickness, and y 0 is one-half the plate spacing. (Swift, 
1988, p. 1157) 

From the equation of motion for a compressible, viscous 
fluid. Swift has shown that 




( 10 ) 



d£i _ -jcjp m ( ) = (Jpj uj* ) 
dx l - £» l - f „ ' 



(ID 



where 



tanh [ (l + j)y 0 /5J 
(1 + j)y 0 /5 u 



( 12 ) 
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and S v = V 2/Lt/pcj is the fluid's viscous penetration depth and 
p is the dynamic viscosity of the fluid (Swift, 1988, p. 
1159) . 

Finally, the wave equation can now be used to solve for 
d<u,>/dx. (Swift, 1988, p. 1157) 



(i + { r 1)f * )p i + 






1 + e. 



u' 



dx 



( 



4 - A 



(P 2 (1 - O) (1 + € s ) 



P/n 

c? Pl 
r dx 



f v d Pl 

dx 



(13) 



= 0 



where 



p m c p 6 K tanh[ (l » j)y 0 /6J 
p s c s 6 s tanh[ (l + j) l/b a ] 



(14) 



c p and c, are the isobar ic specific heats per unit mass for the 
fluid and the solid (i.e. the plate), = / 2K/p m c p w and 6, = 
/ 2K,/p s c,o) are the thermal penetration depths in the fluid and 
the solid, p s is the density of the solid, and K and K, are the 
thermal conductivity of the fluid and the solid. Also, y is 
ratio of the fluid's isobaric to isochoric specific heats, a 
is the fluid's Prandtl number, VT m is the mean-temperature 
gradient in the x-direction, T m is the mean temperature, and 
f K is a function similar to f u . Substituting Eq. (11) into the 
second term of the wave equation yields the following; 
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( 15 ) 



d u x ) 

dx 



= J 



f K ~ f V 



wp m (1 - a) (1 + e s ) 
• w 



<*Pi 
T dx 



Pm*‘ 



n * 4 ^> p , 



Applying these four expressions, Eqs. (9) , (10) , (11) , and 
(15), to Eq. (8) gives the equation we have been working for; 



w = 



— fle | 
2 



P\ 0) 



[j< 



1 + 



Vr 



P»a‘ 



l - (i - o) (l + e s ) kr ; 



- tan (kx) + 1 



( 16 ) 



1 ,2 



(1 + — ) 

+ cos 2 (kx) - j 



1 + e £ 



1 - f. 



sin 2 (kx) ] 



Now with w known, computing the total acoustic power 
simply requires integrating over the volume of the fluid 
within the stack; 

^ stack = / ” dV = / "Sdx. ( 17 ) 



C. TOTAL ACOUSTIC POWER, W 

To compute 1/Q, the total acoustic power W is needed. In 
the previous section, only an expression for the work done by 
the prime mover's stack was derived. The four remaining 



12 



sections of the prime mover also need to enter into the 
equation. The total work done by the prime mover can be 
divided into the sum of the work done by or on each of these 
five individual sections, 



W ^amb end + ^ 'amb hx + ^ stack + ^ 'hot hx + ^ hot end’ 



(18) 



It should be noted that, with the exception of W sUck , these 
terms represent losses. W sUck can represent either loss or gain 
depending on the value of the temperature gradient. 

The equation for w was derived for the stack region but it 
actually applies to each of the remaining four regions. The 
parameters in w simply vary depending on what section of the 
prime mover x happens to be in. In all four remaining regions 
for example, VT m = 0 since there is no temperature gradient 
anywhere except in the stack region. Likewise, 1 = 0 in the 
ambient and hot ends since these regions do not contain any 
plates and the functions f, and f, take on slightly different 
forms as well. 

The total acoustic power is consequently one integral from 
the beginning of the ambient end to the end of the hot end, 

W = f L w(x ) S(x) dx. (19) 

Jx*0 



Within this integral numerous parameters vary with respect to 
x depending upon which region of the prime mover x happens to 
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be in (e.g., perimeter size, plate thickness, plate spacing). 
Also, within the given region of the stack there are several 
parameters which vary as a function of x because of the 
temperature change along the stack (e.g., sound speed, thermal 



parameters are treated as functions of x so that a single 
integration yields the total work done by the prime mover. 

D. STORED ENERGY, E ST 

The final quantity needed to calculate 1/Q is the energy 
stored in the resonator. This is found by integrating the 
time averaged energy density over the volume of the resonator. 



conductivity of the fluid, viscosity, etc.). All such 




( 20 ) 



where the instantaneous energy density is given by 




( 21 ) 



and its time average by 




( 22 ) 
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Again, as in the case of acoustic power, the total stored 
energy can be represented as the sum of the stored energy in 
each of the prime mover's five regions, 

E qT = E qr + E qr + E qr + E qr + E qr . (23) 

hx stack ^‘‘■bothx hot and 



Choosing appropriate expressions for p, and u, s and integrating 
the time averaged energy density over each of the five regions 
produces : 



^ST 



A ? I amb end 

4 P.'S 21 



f t co3 . g . -L)sin*lkx * *)]<* 

JamJbhx 1 + _1_ y o 

Vo 

f [ CQ3 Hkxj_Vl > (1 , J.) sln * (toc 

JaeacJc i + J_ y o 

Vo 

f ( CosUkx ^) t t J_ )sin 2 Ux t 4, )ldx 

J hot hx -i . 1 



0U 



+ X, 



1 

hot endj 



1 + — 



where <p is the phase of the standing wave at the entrance of 
the region. Finally, combining Eqs. (16), (19) and (24) with 
Eq. ( 6 ) yields the desired parameter, 1/Q. 
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III. EXPERIMENTAL SETUP 
A. THE PRIME MOVER TUBE 

The physical construction of the prime mover assembly was 
previously shown in Figure 1. Its internal diameter is 3.82 
cm and its overall internal length is 1.0005 m. A complete 
listing of all important prime mover measurements is found in 
Appendix A. 

1. The Hot Heat Exchanger and the Hot End 

The hot heat exchanger and the hot end tube are both 
constructed of nickel as is the end-cap of the hot end tube. 
This end-cap is welded to the hot end tube and presents a 
rigid boundary for sound waves generated by the prime mover. 
Through the center of this end-cap a small hole has been 
drilled to accommodate a thermocouple probe for sensing the 
hot heat exchanger temperature. To prevent leakage of gas 
through this hole, the thermocouple probe is epoxied in place. 

The hot end tube is 0.63 cm thick and has an internal 
length of 5.50 cm. Its internal diameter, as are all of the 
individual components which make up the prime mover, is 3.82 
cm. The hot heat exchanger assembly consists of 25, 0.045 cm 
thick and 0.762 cm long nickel plates. The spacing between 
each plate is 0.104 cm. To minimize deformation from high 
temperatures, spacers made from 304 stainless steel with a 
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diameter of 0.031 cm are welded to one side of each plate. 
The hot heat exchanger assembly is mounted to the open end of 
the hot end tube. 

2 . The Stack 

The prime mover stack is housed in a 304 stainless 
steel tube of 0.05 cm thickness with a length of 3.50 cm. It 
is comprised of 35, 0.028 cm thick 304 stainless steel plates 
and is the same length as the tube. The plate spacing in the 
stack region is 0.077 cm and, as in the hot heat exchanger, 
spacers are attached to each plate to minimize deformation. 
The stack and the hot heat exchanger are assembled such that 
their plates are at near right angles to one another so that 
the plates of the heat exchanger do not coincide with the gaps 
in the stack. Such a misalignment would result in the 
blockage of fluid flow. 

3. The Ambient Heat Exchanger and Ambient End 

The ambient heat exchanger consists of two identical 
copper stacks each containing 25, 1.02 cm long, 0.045 cm thick 
plates. The gap between each plate is 0.104 cm and the 
spacing between the two individual stacks is 0.15 cm. These 
two stacks are housed within a brass water jacket and the hot 
end/stack assembly is soldered to one end of the jacket. The 
other end of the jacket is soldered to the ambient end of the 
prime mover which is constructed of an 88.11 cm long copper 
tube. A butterfly valve is positioned within this tube, 38 cm 
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from the ambient heat exchanger. The purpose of this valve is 
to shut off the sound wave generated by the prime mover. By 
closing the valve the internal length of the prime mover is 
effectively shortened requiring a much higher differential 
temperature to drive the prime mover above onset. In this 
way, a differential temperature which places the prime mover 
above onset can be established without sound being generated. 
By opening the valve the prime mover is instantaneously placed 
above onset (turned on) and the initial pressure wave build up 
may be observed and measured. Lastly, the ambient end is 
flanged and a brass end-cap which houses a pressure sensing 
transducer is fastened to it. 

B. CONTROL AND MEASUREMENT EQUIPMENT 
l. Temperature 

The temperature of the hot heat exchanger is 
controlled by an Omega Engineering HBA Model 202040 heater 
powered by a GR Type 100-Q variable transformer. The heater 
coils are tightly wrapped and cover the entire hot end and hot 
heat exchanger assembly to provide a constant "hot" 
temperature to this portion of the prime mover. As previously 
mentioned, a type K thermocouple entering through the hot end 
end-cap is positioned in contact with the hot heat exchanger 
to measure the hot end temperature. During operation the hot 
end assembly is surrounded with insulation to minimize heat 
loss . 
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Ambient end temperature is controlled through the use 
of a Neslab RTE-110 constant temperature bath. The bath is 
used to cycle constant temperature water through the water 
jacket surrounding the ambient heat exchanger and then through 
plastic tubing which is wrapped around the entire length of 
the ambient end. Temperature is monitored via three type E 
thermocouples which are glued to the top, middle, and bottom 
sections of the ambient end. This method of ambient 
temperature control maintains the ambient end temperatures 
within a 1°C spread. The hot end thermocouple along with the 
three ambient end thermocouples are monitored using a Keithley 
740 System scanning thermometer. 

2. Pressure 

Mean pressure inside the prime mover is sensed via a 
connection which penetrates the prime mover between the two 
ambient heat exchanger stacks. This connection is valved so 
that the pressure sensing line may be isolated from the prime 
mover when oscillations are being measured. Leaving the valve 
open increases the losses of the prime mover and actually puts 
it below onset for the differential temperatures examined. 
Pressure is measured from this sensing line by an Omega PX304- 
300AV pressure transducer and an HP-3478A multimeter is used 
to display the voltage output of the transducer. Pressurizing 
the prime mover is done by first evacuating the prime mover 
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chamber with a vacuum pump and then filling it with helium 
gas. All this is done via the pressure sensing line. 

3. Signal 

The acoustic signal generated by a differential 
temperature applied to the prime mover is sensed by an ENDEVCO 
Model 8510B-5 piezoresistive pressure transducer housed within 
the ambient end end-cap. The linear response range of the 
transducer is ±5 psi. A high impedance leak is provided 
between the resonator and the back volume of the transducer to 
eliminate the dc pressure difference on the transducer while 
having negligible effect on the acoustic pressure differences. 
The leak is a 0.004” internal diameter, approximately 10 cm 
long piece of copper-nickel capillary tubing. The signal from 
the transducer is routed through a locally constructed 
differential amplifier which boosts the signal by a factor of 
ten. This amplified signal is then applied as input to a 
Nicolet 310 digital oscilloscope for data collection and 
storage . 
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IV 



EXPERIMENTAL PROCEDURE AND SIGNAL PROCESSING 



A. DATA COLLECTION 

The first step in obtaining prime mover response data was 
to evacuate the resonator. This was accomplished by pumping 
down the prime mover via its pressure sensing line for several 
hours. Once fully evacuated, the vacuum pump was turned off 
and isolated from the prime mover. The tube's internal 
pressure was then recorded and this value was used as the 
calibration value for the pressure sensor. All future 
pressure readings were corrected by this amount. The prime 
mover was then filled with the desired pressure of helium to 
within a tolerance of ± 1 kPa. 

For temperature control, the constant temperature bath was 
started and the variable transformer was turned up to power 
the heater element. Once a steady state differential 
temperature of interest was reached, the prime mover was ready 
for data collection. 

As mentioned previously, a Nicolet 310 digital 
oscilloscope (NIC-310) was used to sample the waveform 
generated by the prime mover. Since the data was being 
digitally sampled it was important to at least satisfy the 
Nyquist sampling rate of 2f nux in order to prevent aliasing, 
where f^ is the highest frequency of interest (Strumm, 1988, 
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p. 54) . For the prime mover, we are only concerned with the 
rise time of the fundamental component of the acoustic signal 
which has a frequency of approximately 500 Hz. Consequently 
we must sample at a minimum of 1000 Hz. To be conservative a 
sampling frequency of four times the fundamental was chosen as 
the minimum acceptable sampling frequency. As it turns out, 
the closest sampling frequency offered by the NIC-310 without 
violating this threshold is 5000 Hz or 200 ns per data point. 
With the 4000 point data capacity of the NIC-310, this process 
enabled the recording of sample records of up to 0.8 seconds 
in duration. 

Data collection was done by placing the onset control 
valve in the closed position to turn off the sound generated 
by the tube. After thermal equilibrium was achieved, the 
valve was opened. When a signal could be seen over the 
background noise present, the oscilloscope record was frozen 
and saved on a 3.5" floppy disk for later analysis. Several 
such records were taken for each pressure and temperature 
combination of interest. 

B. DIGITAL SIGNAL PROCESSING 

Analysis of the accumulated data was done with MathWorks' 
program MATLAB. The actual data records stored by the NIC- 
310, however, were in binary code so, to make the data records 
ready for processing, they were simply converted into ASCII 
files . 
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Initially, for each data record, a Discrete Fourier 
Transform ( DFT) was computed and the fundamental frequency of 
the record was identified. This frequency was then used to 
translate the frequency spectrum of the signal so that the 
fundamental frequency was now located in the first bin of the 
DFT. This was accomplished by multiplying the signal by a 
complex exponential in the following fashion; 

x new (n) = e j{2n/lf)Mn x old {n) , (25) 

where N is the number of points contained in the data record, 
n is an integer denoting the position within a record, and M 
is the number of bins which every frequency bin in the DFT of 
x old (n) is to be circularly shifted. For example, if the 
fundamental frequency of the DFT of x old (n) resides in bin M + 
1, then a DFT of x^fn) will have the fundamental frequency 
located in its first bin. (Strumm, 1988, p. 429) The DFT of 
a typical data set is displayed in Figure 3. Note that there 
are two peaks in this DFT. This is because it is the DFT of 
a real (non complex) signal. Since real signals are 
sinusoidal, (i.e. the sum of two exponential functions of 
equal magnitude, but with frequencies of opposite sign) , two 
fundamental frequencies result. One represents the positive 
frequency of the sine wave, the first half of the DFT, and the 
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Frequency Bin 



Figure 3 Discrete Fourier Transform (DFT) of a typical 
data record showing the presence of a 500 Hz signal. 
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other represents the negative frequency of the sine wave, the 
second half of the DFT. 

Now, with the fundamental frequency effectively translated 
to dc level (first bin), a low pass filter can be applied to 
isolate the desired signal. The low pass digital filter 
chosen for this task was a fifth order Butterworth filter with 
a bandwidth of 50 Hz. This procedure is essentially the same 
as filtering the original data with a 100 Hz band pass filter 
centered at the fundamental frequency of the data record. By 
translating the data first, however, we immediately arrive at 
the envelope of the waveform since its frequency has been 
"converted" to dc. Figures 4 and 5 show a typical data record 
before and after processing, respectively. Note that 
processing the signal in this manner has the effect of halving 
the magnitude of the original signal. 

In most cases the 4000 point DFT revealed that only the 
fundamental frequency was present. As the prime mover was 
driven higher (-50K) above onset, however, the second and even 
the third harmonics were observed. To ensure that only the 
affects caused by the fundamental were measured, DFT' s were 
taken of the earliest parts of these data records to eliminate 
the second harmonic. This is important because, as the second 
and higher harmonics begin to grow, they take energy away from 
the first harmonic and cause /3 to change. Such losses are not 
accounted for by our theory. This effect can be seen in 
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Pressure = 500 kPa, Temperature Difference = 388.7 K 




Figure 4 Typical data record before processing 
(Sampled at 5000 samples per second.) 



xlO -3 Pressure = 500 kPa, Temperature Difference = 388.7 K 




Figure 5 The result of processing the data record 
shown in Figure 4 above. 
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Figures 6 and 7 where the first 2000 data points contain only 
the fundamental frequency while the entire 4000 point sample 
shows evidence of the first three harmonics. Note that for 
the first 0.4 seconds (2000 samples), the slope of the 
logarithm of the envelope is essentially constant (Figure 7). 
This slope is equivalent to jS since d/3/dt is approximately 
zero. As the harmonics begin to appear, however, the 
magnitude of /? noticeably decreases as indicated by the 
lessening of the slope on the logarithmic plot of the data. 
The 2000 and 4000 point DFT's of this data are shown in 
Figures 8 and 9. 

As with all digital filtering, the end conditions are the 
most difficult to match. Because of this the initial and 
final points of a filtered signal need to be discarded. To 
determine the actual range of the filtered data points which 
are "influenced" by the filtering process, test signals were 
generated of the following form; 

TestSignal = 0.001 e pt sin(2* 500 t) . < 26 ) 

The use of a scaling factor of 0.001 and a frequency of 500 Hz 
were made to generate signals which were characteristic of the 
data being analyzed. Likewise, varying values of (3 were used 
to cover the range of /?'s that were being measured. The 
results show that smaller values of j8 produce larger relative 
fluctuations in the filtered waveform. Figures 10 and 11 show 
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Figure 6 Filtered signal (solid line) which 
experiences the growth of higher harmonics after 
the first 0.4 seconds. The dashed line is a least 
mean squares fit to the first 0.4 seconds of data. 




Figure 7 Logarithmic plot of signal shown in 
Figure 6 (solid line) . The dashed line is a least 
mean squares fit to the first 0.4 seconds of data. 
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Figure 8 A 2000 point DFT of the signal shown in 
Figures 6 and 7, showing that only the fundamental 
frequency is present in the first 0.4 seconds of 
the signal. 




Figure 9 A 4000 point DFT of the signal shown in 
Figures 6 and 7, showing that three harmonics are 
present in the entire 0.8 second sample. 
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xlCH Test Signal with Beta = 0.5 




Figure 10 A filtered test signal with /3 = 0.5 

showing the endpoint influence caused by the 
filter. 



Test Signal with Beta = 0.5 




Figure 11 A logarithmic plot of Figure 10 showing 
the endpoint influence of the filter on a smooth 
test signal. 
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both the linear and logarithmic presentations for a test 
signal (/? = 0.5) that was processed in accordance with the 
above procedure. As these plots indicate, the filter's 
endpoint effect extends to the first and last several hundred 
data points. To verify this, a least mean squares fit was 
applied to the natural logarithm of the test signal shown. In 
an attempt to eliminate the filter's endpoint effect the first 
and last 300 points were discarded. The result was that (3 
could be determined to within one millionth of one percent, a 
negligible error. For data records which contained only the 
fundamental harmonic, this criteria of discarding the initial 
and final 300 data points, after filtering, was used to 
prevent the filter from adversely affecting the processed 
results . 

In the case where a signal was observed to contain 
multiple harmonics, the filtered data was truncated at some 
value less than 3700. The actual value was determined by two 
criteria. First, as previously mentioned, the DFT was used to 
determine at what point in a data record, the buildup of the 
second harmonic began. This value marks the maximum 
acceptable truncation level. A least mean squares fit was 
then applied to the logarithm of the filtered signal truncated 
at this value. To evaluate the quality of this truncation, 
the filtered signal was plotted along with the least mean 
squares fit conducted over this truncated range. If the two 
curves did not visibly agree over the entire truncated range 
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then it was obvious that the fundamental was not completely 
isolated. This can occur because the initial buildup of a 
harmonic is difficult to detect. Because we are dealing with 
very narrow band frequencies, they will initially go 
undetected since the frequency bins are 1.25 hertz wide for a 
0.8 second sample. Also, as less data is used in performing 
a DFT in an attempt to isolate the fundamental frequency, the 
frequency bins become wider. In cases where the results of 
the initial truncation were unsatisfactory, a smaller portion 
of the sample was chosen and the data was reprocessed under 
this same criteria until an agreement was evident. Appendix 
B gives a tabulation of each data record analyzed and the 
range of the points that were used for determining /3. 

This chosen segment of the filtered data record now 
represents the exponential behavior of the fundamental 
frequency alone and, assuming that d/3/dt is small, the natural 
log of this data is a line whose slope is equivalent to /3. By 
performing a least mean squares fit to this data, a value for 
j3 can be obtained. Figures 12 through 15 show four typical 
data records. They reveal that £ does not noticeably change 
over the time that the data is observed when only the first 
harmonic is present. The uncertainty in the slope was 
determined according to a procedure discussed in Beers' theory 
of errors (Beers, 1957) . The maximum measurable error in all 
cases was only a small fraction of one percent. 
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Pressure = 238 kPa, Temperature Difference = 339 6 K 




Figure 12 A logarithmic plot of a typical data 
record (solid line) . The dashed line is a least 
mean squares fit of the data displayed. 



Pressure = 308 kPa, Temperature Difference = 349.4 K 




Figure 13 A logarithmic plot of a typical data 
record (solid line) . The dashed line is a least 
mean squares fit of the data displayed. 
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Figure 14 A logarithmic plot of a typical data 
record (solid line) . The dashed line is a least 
mean squares fit of the data displayed. 




Figure 15 A logarithmic plot of a typical data 
record (solid line) . The dashed line is a least 
mean squares fit of the data displayed. 
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V. RESULTS AND DISCUSSIONS 



A. ABOVE ONSET DATA 

Data was obtained for a helium filled prime mover above 
onset for mean gas pressures of 238, 308, 376, and 500 kPa , 
and differential temperatures ranging from onset up to 
approximately 400 K, in 10 degree increments. The pressures 
were chosen because they cover the range of the "ideal" gas 
pressure for the plate spacing in the prime mover stack and 
because they match the pressures analyzed in Lin's 1989 
thesis. Although Lin also analyzed the pressure of 170 kPa, 
it was not used in this thesis. This is because the onset 
temperature of the prime mover at 170 kPa is so high that an 
insufficient range of data could be measured with a self 
imposed hot end temperature maximum of 700 K. 

The results for the four pressures are presented in 
Figures 16 through 19 which show 1/Q as a function of the 
temperature difference across the stack. As can be seen, the 
collected data closely follows the slope of the theoretical 
curves and, for high mean gas pressure, the data and theory 
exactly agree. As the mean gas pressure decreases, however, 
the actual value at which onset occurs becomes increasingly 
higher than the theoretical value. Consequently, the spread 
between the actual data and the theory increases as well. 
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Figure 16 Data obtained for 238 kPa plotted as a 
function of the temperature difference across the 
stack. The solid line is the theoretical result. 




Figure 17 Data obtained for 308 kPa plotted as a 
function of the temperature difference across the 
stack. The solid line is the theoretical result. 
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Figure 18 Data obtained for 376 kPa plotted as a 
function of the temperature difference across the 
stack. The solid line is the theoretical result. 




Figure 19 Data obtained for 500 kPa plotted as a 
function of the temperature difference across the 
stack. The solid line is the theoretical result. 
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The results can also be examined from a different point of 
view. Instead of an absolute comparison of theory and data, 
it is useful to consider 1/Q as a function of the number of 
degrees above onset. In this way, the actual results and the 
theory simply start from the common reference point of zero 
quality factor. This will consequently eliminate any 
discrepancies that result from failure to predict the precise 
onset temperature. Figures 20 through 23 present the data and 
the theory in this fashion. To determine the onset 
temperature for the measured data, a least mean squares fit 
was performed on the data and the result was extrapolated to 
find where 1/Q equaled zero. What these plots reveal is that 
at high mean gas pressures, theory exactly agrees with the 
measured data over the entire range of differential 
temperature measured (approximately 80 K) . Again however, as 
the mean gas pressure decreases, the data visibly deviates 
from the expected theory within this differential temperature 
range observed. 

Summarizing, the data displays two effects relating to the 
mean gas pressure. First, as the mean gas pressure decreases, 
the theoretical value of onset is increasing less than the 
measured value; and second, as the mean gas pressure 
decreases, the slope of the actual data is increasingly less 
than the theoretical slope. Three possible explanations exist 
as to why these deviations occur. First, imperfections in the 
construction of the prime mover may cause the prime mover's 
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Figure 20 Data obtained for 238 kPa plotted as a 
function of the number of degrees above onset. The 
solid line is the theoretical result. 




Figure 21 Data obtained for 308 kPa plotted as a 
function of the number of degrees above onset. The 
solid line is the theoretical result. 
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Figure 22 Data obtained for 376 kPa plotted as a 
function of the number of degrees above onset. The 
solid line is the theoretical result. 




Figure 23 Data obtained for 500 kPa plotted as a 
function of the number of degrees above onset. The 
solid line is the theoretical result. 
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output to be less than ideal. Second, a difference in the 
rate of harmonic growth could be responsible for the data 
which is deviate from the theory. And third, the deviate data 
could have resulted if this data was measured during 
conditions of higher thermoacoustic heat transport. 

Since the prime mover is not of ideal construction, its 
performance would be expected to be less than the theoretical 
prediction. This is what has been observed. In all cases, 
the measurements, at best, matched the theoretical results. 
The actual energy transport taking place in the prime mover 
stack occurs in a fluid channel which is approximately one 
thermal penetration depth thick (<SJ and centered 
approximately one thermal penetration depth from the plate 
(Swift, 1988, p.1152). This is important because as the mean 
gas pressure is decreased, 5. increases. Since <S, is also a 
function of temperature, its value changes along the stack. 
For the two lowest pressures looked at though, S K at the 
average stack temperature, is greater than the half-gap 
spacing of the plates, and for the two highest pressures, it 
is smaller. Consequently, any imperfections in the 
arrangement of the plates would increasingly interfere with 
the prime mover's performance as the mean gas pressure is 
decreased. At higher mean gas pressures, the misalignment of 
the plates simply becomes less of a concern because the heat 
transport occurs closer to the plates. 
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The theory does not take into account energy lost to 
harmonics. So if harmonics are present, they would result in 
the absolute value of 1/Q being less that predicted. Analysis 
reveals, however, that harmonic growth is not the cause of the 
deviations. This is evident because the growth of higher 
harmonics was observed to be relatively independent of the 
mean pressure. Also, harmonic generation should be stronger 
for larger acoustic amplitudes. Therefore, it would 
reasonably be expected that the higher AT data would be more 
effected than the lower AT data since /3 is much larger. It 
can be seen however, from Figures 20 and 21, that the slope of 
the measured data, though in deviation with the theory, is 
very linear. No AT dependence is evident. Also, we were 
careful to exclude data records which had noticeable harmonic 
content. 

Lastly, the question of heat transport must be addressed. 
The entire reason for looking at signals of low amplitude is 
to ensure that the heat transported by the sound is not yet 
significant. This is desired because, after the heat 
transport becomes significant, the heat exchangers initially 
have trouble keeping up with the energy demand placed upon 
them. The transport of heat through the fluid is simply 
faster than the transfer of energy through the solid heat 
exchangers. Consequently, the temperature of the hot and 
ambient end begin to change and the stack's condition becomes 
dynamic. With thermal equilibrium lost, the measured values 
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no longer represent / 3 . Since heat transport is proportional 
to the square of the amplitude of the pressure wave, it is 
important to know at what pressure the signal was at when 0 
was measured. Figures 24 through 27 show the range of the 
peak acoustic pressures generated by the prime mover in the 
data records used in determining the values of / 3 . The 
pressures are small enough that the resultant heat transport 
is insignificant. Also, they are all approximately the same 
order of magnitude for the various mean pressures. As a 
result, any heat transport would be essentially the same for 
high and low mean pressures. Finally, deviations due to heat 
transport would be AT dependent because the pressure 
amplitudes are higher for higher values of AT. As noted 
above, however, no AT dependence is evident. It can be 
concluded therefore, that heat transport is not responsible 
for the deviations. 

B. SOURCES OF ERROR 

The theory discussed in Chapter II does not easily handle 
the insertion of the butterfly control valve in the prime 
mover ambient end. To see how much effect the valve has on Q, 
the Q of the prime mover was measured with no temperature 
difference applied to the stack. These results were then 
compared to Lin's measurements of Q for the same prime mover 
but without the valve. The measurement procedure used to 
obtain Q was identical to the procedure used by Lin in 1989. 
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Mean Pressure = 238 kPa 




Figure 24 Range of the peak acoustic pressure of 
each data record used in determining /3. (P = 238 
kPa) 



Mean Pressure = 308 kPa 




Figure 25 Range of the peak acoustic pressure for 
each data record used in determining /3. (P = 3 08 
kPa) 
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Mean Pressure = 376 kPa 




Figure 26 Range of the peak acoustic pressures for 
each data record used in determining 0. (P = 376 
kPa) 



Mean Pressure = 500 kPa 




Data Record 

Figure 27 Range of the peak acoustic pressures for 
each data record used in determining /3. (P = 500 
kPa) 
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The results of this comparison show that the Q's agree within 
-1%. (Lin, 1989, pp. 21-34) 

The largest experimental error is likely due to the 
measurement in temperature. The error, based on 

manufacturer's specifications, is ±2°C for type E 
thermocouples and ±5°C for type K over the temperature range 
of interest. 

As discussed at the beginning of Chapter IV, the error in 
mean pressure is ±1 kPa. 
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VI 



CONCLUSIONS 



It has been observed that the theoretical predictions of 
the work output of a prime mover are in very good agreement 
with the data for high mean gas pressures. As mean gas 
pressure decreases, however, the agreement diminishes. 
Likewise, the prediction of the slope of the data above onset 
is very good at high mean gas pressures. It is exact over 
large ranges of differential temperature. Again however, as 
the mean gas pressure decreases, the slope of the theoretical 
values of 1/Q increasingly exceed those of the experimental 
results. These results are consistent with those found for 
the same prime mover below onset (Lin, 1989) . As explained in 
the previous chapter, these observed deviations are probably 
due to the inevitable construction irregularities of the prime 
mover stack. This method for measuring 1/Q has proven to be 
very effective for low amplitude signals. 

To verify the discussion of the results in the previous 
chapter, it is recommended that one of two things be done. 
First, a working gas other than helium could be chosen. As a 
result, the thermal penetration depth S K would be different 
from those employed in this thesis. If the same S K to plate 
spacing ratios are examined as those in this thesis, then 
similar results are expected. Second, a new prime mover stack 
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could be constructed with a plate spacing different from that 
of the stack used here. When analyzed at pressures that give 
the same penetration depth to plate spacing ratio, similar 
results are expected. If either experiment produces results 
that differ from those seen here, then it would indicate that 
the source of the deviations observed in this thesis is 
something more than just a physical construction problem. 
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APPENDIX A 



PRIME MOVER DIMENSIONS 



Prime Mover Internal Length 


Prime Mover Internal Radius 


100.05 cm 


1.91 cm 



COMPONENT DIMENSIONS AND LOCATIONS 



Component 


Ambient Heat 
Exchanger 


Stack 


Hot Heat 
Exchanger 


Width (Ax) 


2 . 18* cm 


3.50 cm 


0.762 cm 


Plate half- 
gap (y 0 ) 


0.052 cm 


0.0385 cm 


0.052 cm 


Plate half- 
thickness (1) 


0.0225 cm 


0.014 cm 


0.0225 cm 


Total plate 
perimeter 


3 12 cm 


206 cm 


156 cm 


Location (x) 


88.11 cm 


90.29 cm 


93.79 cm 



* Includes 0.15 cm gap between twin heat exchangers. 
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APPENDIX B 



DATA RECORD SUMMARY 



, 

! Pressure (kPa) 


Temperature 
Difference (K) 


Range 

containing 
only first 
harmonic 


Filtered 
range used 
in jS 

measurement 


238 


339.6 


1-4000 


300-3700 


238 


349.2 


1-3000 


300-3000 


238 


359.7 


1-4000 


300-3700 


238 


367.2 


1-4000 


300-3700 


238 


379.6 


1-3000 


300-2900 


238 


388.1 


1-3000 


300-3000 


238 


349.6 


1-2000 


300-2000 










308 


332.5 


1-4000 


300-3700 


308 


342.7 


1-4000 


300-3700 


308 


349.4 


1-4000 


300-3700 


308 


357.6 


1-4000 


300-3900 


308 


366.3 


1-4000 


300-3700 


308 


378.4 


1-4000 


300-3700 


308 


369.4 


1-3900 


300-3000 


308 


397.0 


1-2000 


300-1500 
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DATA RECORD SUMMARY (CONTINUED) 



Pressure (kPa) 


Temperature 
Difference (K) 


Range 

containing 
only first 
harmonic 


Filtered 
range used 
in /? 

measurement 


376 


339.3 


1-4000 


300-3700 


376 


360.3 


1-3000 


300-3700 


376 


360.3 


1-3000 


300-3700 


376 


367.8 


1-3000 


300-3700 


376 


378 . 9 


1-3000 


300-3500 


376 


387.9 


1-4000 


300-2000 


376 


398 . 1 


1-3000 


300-2000 










500 


361.7 


1-3000 


300-3700 


500 


37C.9 


1-4000 


300-3700 


500 


379 . 6 


1-4000 


300-3700 


500 


388.7 


1-3000 


300-3700 


500 


398 . 2 


1-4000 


300-3700 | 
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APPENDIX C 



EQUATIONS USED FOR DETERMINING 
TEMPERATURE DEPENDENT QUANTITIES 

Mean Temperature as a function of x : 



k T0 


+ 


krpi ( x x sCack ) + k T2 (x x sCack 


k T0 


= 


C T00 + C T01^T 


k T i 


= 


C T10 + + C T12 (AD 2 


k T2 


= 


T20 ^ ^ T22 ^ ^ 


C T00 


= 


292.86 


C T01 


= 


0 . 00171124 


C T1 0 


= 


-12.6165 


C T11 


= 


29.0362 


C T12 


= 


0.00920377 


C T20 


= 


335.803 


C T21 


= 


-12.3295 


C T22 


= 


-0.270547 



Mean Temperature Gradient : 

Vr ro [x] = k g0 + k gl e~ k ° lU ~ <*' 



kg() 


= 


c goo + c goi^T + 


^g02 ( A T) 


k 9l 


— 


CglO + C gll AT + 


C gl2 ^ ^ ^ 


k g2 


= 


c g20 + C g21^T + 


^ g2 2 ( ^ T) 


C g00 


= 


9.36469 




C g01 


= 


18.9501 




Cg 02 


= 


0.0000435077 




C glo 


= 


-9.90268 




c gli 


= 


9.63328 




C gi2 


= 


0 . 0122971 





c g20 = 0 • 121556 
c g21 = 0.07 2249 2 
c g22 = _0 • 0000210737 



Sound Speed: 

a[T] = a 0 JWT 0 

a 0 = 1008 m/ s 
T 0 = 29 3 K 



Angular Frequency : 

»[An = 



rm 



^3 tack 



W 1 IXscack+^XsCacJ 



(1 - X, 



stack 

atzwi 



Ax 



stack 1 



«1 [x] - * vrix] u - X„, ck ) - 



Phase : 

<Mx] = ^[AT]^ " vrixj (x - x sCacJt ) - /T^] 
<MAr] - 2w[An/T^/(a 0 Vr[x] ) 

Miscellaneous : 

/c[r] = o) [AT] /a [T] 

H[T] = \i 0 (T/T 0 )* 
k[T] =k 0 (T/T 0 )P 
p m [T] = p 0 {P/P 0 ) (r o /D 
5 K [T] = ^2k[t] / (p ”1 rj c Po co i a rj ) 

6 v [r] = ^lilTj/tp^LTJwlATJ ) 
o[T] = (6 V [T] /6 k [T] ) 2 
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Miscellaneous (continued) 



P 

Y 

Ho 

K o 



P O 
Po 



= 0.694 
= 5/3 

= 1.961xl0' 6 Pa-sec 

= 0.1524 W/irrK 

= 0 . 164 kg/m 3 

= 1 . 01x10 s Pa 

= 5197 J/kg-K 
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